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Powers of the space forms curvature operator and 
geodesics of the tangent bundle.* 

Saharova Y., Yampolsky A. 


Abstract 

It is well-known that if T is a geodesic line of the tangent (sphere) 
bundle with Sasaki metric of a locally symmetric Riemannian manifold 
then the projected curve 7 = tt o F has all its geodesic curvatures 
constant. In this paper we consider the case of tangent (sphere) bundle 
over the real, complex and quaternionic space form and give a unified 
proof of the following property: all geodesic curvatures of projected 
curve are zero starting from ^ 3 , fee and fcio for the real, complex and 
quaternionic space formes respectively. 
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Introduction 

Sato K. [4] and Sasaki S. [3] proved that the projection to the base space 
of any non-vertical geodesic line on the tangent or the tangent sphere bun¬ 
dle of a real space form M^(c) is a curve of constant curvatures ki and k 2 
and zero curvatures k^,..., kn-i- Nagy P. [2] essentially generalized this re¬ 
sult. He considered the case of general locally symmetric base manifold and 
have proved that the geodesic curvatures of projection of any (non-vertical) 
geodesic line on the tangent sphere bundle are all constant. Nevertheless it 
was still interesting to find a clearer description of projections of geodesics 
for the case of classical rank one symmetric spaces. The second author made 
a first step in this defection and proved that the projection to the base space 
of any non-vertical geodesic line on the tangent or tangent sphere bundle of 
a complex space form CP^ is a curve of constant curvatures ki,... ,k 5 and 
zero curvatures k^,..., kn-i 

In is this paper we make a contribution in more clear understanding of 
geometry of projected geodesics in the case of tangent (sphere) bundle of 
almost all classical locally symmetric spaces, namely spheres, complex and 
quaternionic projective spaces and their non-compact dual from a unified 
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viewpoint using the recurrent properties of powers of the curvature operator 
of these spaces. This approach allows to give also a unified proof of the 
results from [3], [4] and [5] 

We also use an easy to prove result [1], stating that the geodesics of tan¬ 
gent or tangent sphere bundle with Sasaki metric have the same projections 
to the base manifold. 

Remark on notations. Throughout the paper (•, •) and | • | mean the 
scalar product and the norm of vectors with respect to the corresponding 
metrics. 

1 Summary of main results. 

Let (M"(c),(/) be a Riemannian manifold of constant curvature c, 

(M^"(c); J;g) a Riemannian manifold with complex structure J of constant 
holomorphic curvature c and (M^" (c); Ji, J 2 , J 3 ; g) a Riemannian manifold 
with quaternionic structure {Ji, J 2 -, Jz) of constant quaternionic curvature 
c. For the sake of brevity, denote by M{c) one of these space forms with 
corresponding standard metrics and will refer to A4(c) just to a space form 
of constant curvature c. The main result is the following statement. 

Theorem 1.1 Let Ai{c) be a space form of constant curvature c / 0. Let 
r be non-vertical geodesic line on the tangent or tangent sphere bundle over 
Ai{c). Let 7 = TT o r be the projection of T to A4(c). Then the geodesic 
curvatures ki,k 2 , ■ ■ ■ of ^ are all constant and 

(a) ks = ... = kn-i = 0 for the real space form; 

(b) ke = ... = k 2 n-i = 0 for the complex space form; 

(c) fcio = ... = /c 4 „_i = 0 for the quaternionic space form. 

As the referee remarked, the result of the Theorem 1.1 can be expressed 
in more clear geometrical terms, namely the projected curve 7 = tt o T lies 
in a totally geodesic or , in a totally geodesic CP^ or CH^ and in a 
totally geodesic QP^ or QH^ for the real, complex and quaternionic space 
form respectively. These assertions can be derived from ( 6 ), (10) and (14). 

Proof of the Theorem 1.1 is based on the recurrent property of powers of 
curvature operator of spaces under consideration. Let Rxy be the curvature 
operator of M(c). Define a power of curvature operator R^y recurrently in 
the following way: 


RP^yZ = RP-^{RxyZ) p>l. 

The basic tool for our considerations are a chain of lemmas. 
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Lemma 1.1 Let Rxy be the curvature operator of the real space form 
{M^{c),g). Then for any X and Y 


RP - 

Hxy — 


s > 1 


{—rc^Y ^Rxy for p=2s-l 
{—b"^c^Y~^R \y for p=2s, 

where b = |X A y| is a norm of bivector X AY. 


Lemma 1.2 Let Rxy be the curvature operator of the non-flat complex 
spaee form {M^{c)] J; g). Denote by b = |X A y| the norm of a bivector 
X AY and m = {X, JY). Then for any X and Y 


RYxy — 


Lin {JR?xyi Rxy, J ) for p=2s-l 
Lin JRxy, E) for p=2s, 


s> 2 


where E is the identity operator and Lin means a linear eombination of 
eorresponding operators with coefficients being polynomials in ^ ,b ,m. 

Lemma 1.3 Let Rxy be the eurvature operator of the non-flat quaternionic 
space form {M^{c)] Ji, J 2 , Js] g). Denote by b = \X A Y\ the norm of a 
bivector X AY. Set mi = {X,JiY), m 2 = {X,J 2 Y), m 3 = {X^JYT), 
m? = mf + ml + m 3 , J = mi Ji + m 2 J 2 + rn^J^. Then for any X and Y 


RP - 

Hxy — 


Lin { JR\y, Exy^ JE\y^ Rxy,J) for p=2s-l 
Lin (i?^y, LLR^xy, E?xy, J^Exy^E) for p=2s, 


s > 3 


where E is the identity operator and Lin means a linear eombination of 
eorresponding operators with coefficients being polynomials in ^ ,b ,m. 


2 Necessary facts and proof of the main result. 

Let {M'^,g) be a Riemannian manifold and TM^ be its tangent bundle. 
Denote by {u^,... ,u^) a local coordinate system on M”. Then in each 
tangent space of the natural coordinate frame {d/du ^.,..., d/du'^} form 
a local basis. Let ^ be any tangent vector over the given local chart. Then 
f can be decomposed as 

f = f^d/du^ + ...+fYd/du^. 

The parameters ..., ...,form the so called natural induced 

eoordinate system in TM'^. The Sasaki metric line element da"^ with respect 
to this coordinate system is 

= ds"^ + ( 1 ) 
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where ds'^ is a line element of is the covariant differential of ^ with 

respect to Levi-Civita connection on M” and | • | means the norm with 
respect to Riemannian metric on M^. 

The tangent sphere bundle TiM” can be considered as a hypersurface 
in the tangent bundle defined by the condition | ^ | = 1. We will consider 
TiM^ as a submanifold in TM^ with the induced metric. 

With respect to the natural coordinate system, each curve T on TM"^ can 
be represented as r(cj) = |u^((t) ..., ... , with respect 

to the arc-length parameter a and can be interpreted as the vector held 
^(cj) = {a)d/du^ + ■ ■ ■ +C^{a)d/du'^ along the projected curve 7 = ttoT = 
(m^(cj), ..., u’^(cj)). If ^ is a unit vector held then T lies in TiM" and 
represents an arbitrary curve in TiM'^. 

Denote by {') the covariant derivative along 7 with respect to parameter 
a. Then T is a geodesic line on TM"' or TiM” if 7 and ^ satisfy respectively 
the system of equations 


TM^ : 


i' = 

= 0 ; 


TiM” : 


i' = Ri’ii ■, 
e" = -P^i. 


( 2 ) 


where | P and is the curvature operator of based on bivector 

e A 

From (2) it follows that p = const in both cases. Denote by s the arc- 
length parameter on 7 . Then from (1) it follows that 

^ (3) 

SO that 0 < p < 1. According to the latter inequality, the set of geodesics of 
TM^ and TiM^ can be splitted naturally into 3 classes, namely, 

• horizontal geodesics (p = 0 ) generated by parallel (unit) vector helds 
along the geodesics on the base manifold; 

• vertical geodesics (p = 1 ) represented by geodesics on a hxed hber; 

• umbilical geodesics corresponding to 0 < p < 1 . 

In what follows we will consider the properties of projections of umbilical 
geodesics. 

Lemma 2.1 (cf. [2]) Let {M^,g) be a locally symmetric Riemannian man¬ 
ifold and Rxy its curvature operator. Let 7 = vr o T be a projection of 
geodesic line on TM'^ or TiM^ to the base space. Then for the derivatives 
of 7 of order p we have 


7(p) = 

and as a consequence all the geodesic curvatures of j are constant. 
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Proof. The equalities follow from parallelism of curvature tensor of M" 
and the equations (2). Moreover, from the evident identity 

=0 

for all p > 1 , we conclude that = const for all p > 1 and therefore, by 
induction, all the geodesic curvatures of 7 are constant. 

■ 

Proof of Theorem 1.1. Case (a). Denote by ei,..., e^-i the Frenet 
frame of 7 . Using the Frenet formulas for the curve with constant geodesic 
curvatures and keeping in mind (3), it is easy to see that 


^i 2 s l/ 2 /j^/j 2 ...A; 2 ^_ 2 e 2 s-l + Lin {ei, 63,..., 625-3}, 

^{ 2 s ) = (1 _ p‘^ykik 2 ... k 2 s-i 625 + Lin {62,64, ..., 625-2} 


(4) 


for all s > 1 ( with formal setting fco = 1). Setting s = 1,2 in even deriva¬ 
tives, we see that 


7 ( 2 ) = (1 - p'^)ki 62 

-yU) = (1 - p^)'^kik2kz 64 + Lin (62). 


(5) 


On the other hand, applying Lemma 2.1, Lemma 1.1 and Lemma 2.1 again, 
we get 

7 ^^^ = = —b^c^ ( 6 ) 

Substitution from (5) gives 


(1 - p ) kik 2 kz 64 -h Lin (62) = 0 

and therefore /C 3 = 0 , which completes the proof. 

Remark, that 6 ^ is constant along 7 since 

( 62 )' = = 2p^{i\y)-2{y,i)p^ = t). 

Case (b). Denote by ei, ..., 62 n-i the Frenet frame of 7 . Similar to the 
case (a) considerations, the Frenet formulas give 


^(25 1)=(i_^2^5 1/2/j^/j 2 ...A; 25_2 625-1 + Lin {61,63,..., 625-3}, 

^( 2 ^) = (1 - p^ykik2 ... k 2 s-i 625 + Lin {62,64,..., 625-2} 

for all s > 1 . Setting s = 1 , 2 , 3,4 in odd derivatives, we get 

y =(l-p2)V2g^^ 

-)/( 3 ) = (1 - p^y/‘^kik2 63 + Lin (61), 

7(5) = (1 - p'^y/'^ki ... A:4 65 -h Lin (61,63), 

7(7) = (1 _p2y/2^^... fee 67-h Lin (61,63,65). 


(7) 


( 8 ) 
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On the other hand, applying Lemma 2.1, Lemma 1.2 and Lemma 2.1 again, 
we get 

I 7^'^ = = Lm(E^^,^,JE^,^,E)y = Lin (j(^), ,y}, 

Excluding from (9), we come to the equation 

7 ^'^^ = Lin ( 7 *-^^ 7 ^^\ 7 ^ (10) 

Substitution from ( 8 ) imply 

(1 - p^y/‘^ki ...fee 67 +Lin ( 61 , 63 , 65 ) = 0 


and we conclude that fee = 0 which completes the proof. 

Remark, that the coefficients of all linear combinations are constants. 
Indeed, by Lemma 1.2 the coefficients are polynomials in 1/c, 6 = 1^ A 
and m = , J^). The value h is constant along 7 by the same reasons as in 

case (a). The value m is constant along 7 since 

m' = {yjy' = {e,jy + {yjy)^Q. 


Case (c). Denote by ei,..., 64^-1 the Frenet frame of 7. As above, the 
Frenet formulas give 

^(2s-i) = (1 _ p‘^Y-^/‘^kik2 ... k 2 s -2 62s-i + Lin {61,63,..., 62S-3}, 

7(2^) = (1 - p^Ykik 2 ■ ■ ■ k 2 s-i 62s + Lin {62,64,..., 625-2} 

for all s > 1 . Setting s = 1 , 2 , 3 , 4 , 5 , 6 in odd derivatives, we get 
7' = (1-p2)l/2g^^ 

7(3) = (1 - p^Y/‘^kik2 63 + Lin (ei), 

7(5) = (1 - p^f/‘^ki ... A:4 65 + Lin (ei, 63), 

7^^ = (1 - p^y/'^ki ... fee 67 + Lin (ei, 63, 65), 

7(9) = (1 - p^y/'^ki ... fcs eg + Lin (ei, 63, 65, 67), 

7(11) = (1 - p2)ii/2/j^ ... fcio 611 + Lin (ei, 63,65,67, eg). 

Applying again Lemma 2 . 1 , Lemma 1.3 and then Lemma 2 . 1 , we get 
' 7(7) = E\,^i = Lin{E\,^,JE\,^,E\,^,JE^,^,E)y = 

Lin (7(5), J7(4)^7(3)^ 

7(9) = E\,^y = Lin{E\,^,JE\,^,E\,^,JEi:,^,E)y = 

Lin( 7 ( 5 ), ^7(41,7(3)^ j 7 ( 2 )^y)^ 

7(11) = Rio^y = Lin{E\,^,JE\,^,E\,^,JE^,^,E)y = 

Lin( 7 ( 5 ), ^7(41,7(3)^ j 7 ( 2 )^y)_ 
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(14) 


Excluding and from (13), we come to the equation 

^( 11 ) = Lin(7^®\7^'^\7^^^,7^^\7'). 

Substitution from (12) imply 

(1 - ... fcio eii + Lin (ei, 63 , 65 , 67 , eg) = 0 

and we conclude that fcio = 0 which completes the proof. 

Remark, that the coefficients of all linear combinations are constants. 
Indeed, by Lemma 1.3 the coefficients are polynomials in 1/c, 6 = A 

and m = + ^2 + value b is constant along 7 by the same 

reasons as in case (a). The values mi, m 2 , m 3 are all constant along 7 since 

m' = ie, MY = M) + (e', M') - 0 

for i = 1,2, 3. 

3 Proofs of basic Lemmas 

Proof of Lemma 1.1. The curvature operator Rxy of the real space form 
{M^{c),g) has the following expression 

RxyZ = c [(y, Z)X - {X, Z)Y . 

Then 

R\yZ = c [(y, RxyZ)x - (y, RxyZ)y] = [(y, (y, z)x - {x, z)y)x - 

(y, (y, z)x - (y, y)y)y] = [(y, z)(y, y)x - (x, z) |y- 
(y,z)|x|2y +(y,y)(y,y)y] = c2[(y,z)((x,y)x - jypy) + 

(x, z) ((X,y)y - |y| 2 y)] = c [(y, z)RxyX + (y, z)RyxY . 

Therefore, 

R\yZ = c\bY,RxYZ)RxYX + {X,RxYZ)RYxy\ = [((y, Z)(X, y) - 

(X, Z) |y|2) ((X,y)x - \x\^y) + ((y, z)\x\^ - (x,z)(x,y)) ((x,y)y - 
|y| 2 x)] = c3[- (y,z)x(|x| 2 |y |2 _ (x,y)^) + (x,y)y(|x| 2 |y |2 _ 
(x,y) 2 )] =-cH‘^RxyZ, 

where, evidently, 6 ^ = |X|2|y|2 — l^X,Y)‘^ is the square norm of X A y. 
Now we can find the other powers for Rxy inductively. 

■ 

Proof of Lemma 1.2 

The curvature operator Rxy of the complex space form J; g) 

has the following expression 

R^yZ = I [(y, z)x - (X, z)y + (jy, z)jx - {jx, z)jy + 2 (x, jy)jz . 
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Introduce the unit sphere type operator S acting as 


S{Z) SxyZ = {Y, Z)X - (X, Z)Y, 


and the operator S{Z) acting as 


S{Z) SjxjyZ = {JY, Z)JX - {JX, Z)JY. 

Finally, if we denote m = (X, JY), then the curvature operator under con¬ 
sideration takes the form 




S + S + 2m J 


Z. 


(15) 


Since |X A y| = \{JX) A {JY)\, the operators S and S satisfy 

53 = -h'^S, = -b^S, 


where 6^ = |X A Y\‘^. 

In what follows, we need a ’’table of products” for the operators S and 
S. Namely, 



5 

S 

J 

5 

52 

mJS 

JS 


mJS 

52 

JS 

J 

JS 

JS 

-E 


Indeed, 


(16) 


{SS){Z) = SxY SjxjyZ = Sxy[{JY,Z)JX - {JX,Z)JY\ = 

(y, {JY, Z)JX - {JX, Z)JY^X - i^X, {JY, Z)JX - {JX, Z)JY^Y = 
(y, JX) {JY, Z)X + {JX, Z) (X, JY)Y = 

mj\{JY,Z)JX - {JX,Z)Jy'\ = {mJS){Z), 
{SS){Z) = SjxjySxyZ = Sjxjy[{Y,Z)X - {X,Z)y] = 

(^JY, (y, z)x - {X, z)Y^jx - i^jx, (y, z)x - (x, z)y^jy = 
(jy,x)(y,z)jx + (X, z){jx,y)jy = 

mj[(y, Z)X - {X, Z)y] = {mJS){Z), 
{SJ){Z) = SxyJZ = (y, JZ)X - (X, JZ)Y = 

J[{JY,Z)JX - {JX,Z)JY] = {JS){Z), 

{SJ){Z) = SjxjyJZ = {JY, JZ)JX - {JX, JZ)JY = 

J[{Y,Z)X-{X,Z)Y] = {JS){Z), 


and the other entries of the table can be found in a similar way. 




From (16) we see that J{S + S) = {S + S)J and we have 


(5 + Sf = + SS + SS = + mJ{S + S) 

{S + Sf= {S + S)[S^ + + mJ{S + S)] = 

S^ + S^ + SS^ + 552 + mJ{S + 5)2 = 

-525 - b'^S + (55)5 + (55)5 + mJ(5 + 5)2 = 
-62(5 + 5) + mJ(52 + 52 ) + mJ(5 + 5)2 = 
-62(5 + 5) + mJ[(5 + 5)2 - mJ{S + 5)] + 
mJ(5 + 5)2 = (m2 - 62)(5 + 5) + 2mJ(5 + 5)2. 

Thus, 

(5 + 5)3 = Lin (5 + 5, J(5 +5)2) 

On the other hand, setting for brevity Rxy = R, from (15) we derive 


S + S=^R- 2mJ = Lin {R, J), 
(5 + 5)2 = Lin i^R^^jR^E). 


Comparing (17) and (18) we conclude 


(5 + 5)3 


Lin 


Lin (i?, J), JLin {R?, JR, E) 


Un{JR^,R,J). 


On the other hand, from (18)i 


(5 + 5)3 



i?3 + Lin (Ji?2,J). 


So, finally 

i?3 = Un{JR^,R,J). 

It is easy to trace that the coefficients of all linear combinations are 
polynomials in i, 6 , m. To complete the proof we should remark that 


Lin JLin (Ji?2^ J)^ 


R\JR 


= nm ,, 


which allows to find all powers of R inductively. 


Proof of Lemma 1.3 

The curvature operator Rxy of the quaternionic space form 
(M^"'(c); Ji, J2, J3; g) has the following expression 


RxyZ = 


{Y, Z)X - {X, Z)Y + (JiT, Z)JiX - (JiX, Z)JiY+ 

(J 2 T, Z)J 2 X - {J 2 X, Z)J 2 Y + (J 3 T, Z) J 3 X - (J 3 X, Z) J 3 T+ 
2 {X, JiY)JiZ + 2 {X, J 2 T) J 2 Z + 2 (X, J 3 T) J 3 ZI. 
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where Ji, J 2 , J 3 are operators of quaternionic structure 

J1J2 = J3, J2J3 = Ji, J3J1 = J2, Jf = -E, {X, JiY) = -{J^X,Y), i = M 


Introduce the unit sphere type operator S acting as 

S{Z) SxyZ = {Y, Z)X - (X, Z)Y, 
the operators Si{Z) acting as 

Si{Z) Sj,xj,yZ = {J^Y, Z)JiX - {J,X, Z)JiY, i = M, 
and the operator S{Z) acting as 

S{Z) = Si{Z) + S2{Z) + S^{Z). 

Finally, denote mi = (X, JiY), i = 1,3, rn? = mf + + m 3 , J = mi Ji + 

m 2 J 2 + m 3 J 3 . Then the curvature operator under consideration takes the 
form 


RxyZ = 


4L 


S + S + 2J 


Z. 


(19) 


Since |X A y| = |(JiX) A {JiY)\ (i = 1,3), the operators S and Si satisfy 

= -h^S, Sf = -b^Si (i = M), 

where 6 ^ = |X A 

The table of products for the operators S and S is the following one. 



s 

Si 

S2 

S3 

Jl 

J2 

J3 

s 

S'^ 

miJiSi 

1712 J2S2 

msJsSs 

JiSi 

J2S2 

J3S3 

Si 

mi Ji5 

Si 

-mzJzS2 

-m2J2*S'3 

JiS 

J2S3 

J3S2 

S2 

m2J2S 

-maJzSi 

Si 

-miJi^s 

J1S3 

J2S 

J3S1 

S3 

m^JsS 

-m2J2*S'i 

-mi Ji52 


J1S2 

J2S1 

J3S 

Jl 

SiJi 

SJi 

S3J1 


-E 

J3 

-J2 

J2 

S2J2 

S3J2 

SJ2 

S1J2 

-J3 

-E 

Jl 

J3 

S3J3 

S2J3 

S1J3 

SJ3 

J2 

-Jl 

-E 


The expressions for products SSi, SiS, SJi one can find similar to the table 
(16) making formal replacements S ^ Si and J ^ Ji. As concerns the other 
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entries, we have 


(5i52)(Z) = 

SjiX JiySj^X J2yZ = Sjj^xJiY {J2y, Z)J2X — {J 2 X, Z)J2Y 
'jiY, {J 2 Y, Z)J 2 X - {J 2 X, Z)J 2 Y)JxX- 
J JiX, (JsY, Z)J 2 X - {J 2 X, Z)J 2 Y)JiY = 

{JiY, J2X){J2Y, Z),hX + (JiX, J2Y){J2X, Z)JiY = 

Ji [(X, J3Y){J2Y, Z)X - (X, J3X)(J2X, Z)Y = 

-J1J2 \m^{J2Y, Z) J2X - m3(J2X, Z)J2Y] = {-m^J^S2){Z), 


(5i Ji)(Z) = Sj.xj.yJiZ = {JiY, JiZ) JiX - (JiX, JiZ) JiX = 
Ji[(y,z)x-(x,z)y] = (Ji5)(z), 

{SiJ 2 ){Z) = Sj,xj,yJ 2 Z = {JiY, J 2 Z) JiX - (JiX, J 2 Z) JiX = 

Ji[(j 3 X,z)x-(j 3 X,z)y] = 

- JiJ3[(J3y,Z)J3y - {J3X,Z)J3Y] = iJ2S3)iZ) 

and so on. 

From (20) we see that 

(5 + S)J = (S + Si + S 2 + S3){miJi + m 2 J 2 + rnsJs) = 

miJiSi + m2 J2-S'2 + m3 J353 + mi Ji 5 + m2 J 25'3 + m3 J352 + mi Ji 53 + 

m 2 J2*S' + m3J3Si + mi Ji 52 + m 2 J25'i + m3J3S = 

(mi Ji + m 2 J 2 + m3J3)(5 + ^i + 52 + 53 ) = ^^(5 + 5). 


Therefore, the operators (5 + 5) commute and hence for the operator 
c 

i? = -{(5 + 5) + 2 J'} the usual formula for powers can be applied; 

The powers for J' can be found trivially, since 

= mf Jf + mim2(Ji J2 + J2J1) + mim3(Ji J3 + J3J1) + mlJly 
m2m3(J2 J3 + J3J2) + = —mf£' — m^i? — rngif = —m?E, 

where mf + m^ + m 3 = m^. 

As concerns the powers of 5 + 5, the following proposition gives the 
answer. 


Proposition 3.1 The operator 5 + 5 possesses the recurrent property 

(S + 5)® = -2 (6^ + m^){S + Sf - {h^ - m^){S + 5), 

where 6^ = |X A Xp and rn^ = ml + ml + rnl = (X, JiX)^ + (X, + 2 ^)^ + 
{X,J^Yf. 
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Proof. The proof is technical and in what follows we will use some auxiliary 
operator products. Namely, 


SS = SJ, SS = JS, 

SJS = -m'^S, SSJ = -m'^S, 

s{sf + + Si) = s^j, 552 = JS^ 

SSJ = JSJ, SJS = JS^, ^ ’ 

52 = ,Si 2 + si + si-sj + JS, 

s{sf + si + si) = -b‘^s - {sf + si + si)j + js^. 

The proof is straightforward. Applying (20), we get 
SS = S{Si + 52 + 53 ) = mi Ji5i + m 2 52^2 + m 3 = 
miSJi + m2SJ2 + rrisSJs = SJ. 

In a similar way we find 

SS = {Si + 52 + 53)5 = mi Ji5 + m 2 J 25 + m^J^S = JS, 

52 = (5i + 52 + 53 ) (5i + 52 + 53 ) = 52 + 52 + Si + S 1 S 2 + 5i 53 + 525i + 
5253 + 5351 + 5352 = Sl + Si + Si-m^SiJii-m2SiJ2-miiS2h-miS2J2- 
m2S^J2 - miS^Ji = 5| + 52 + 52 - SJ + mi5iJi + m252J 2 + = 

Sf + 51 + 532 -SJ + JS, 

SJS = 5(miJi + m 2 J 2 + m 3 J 3 )S = (rriiJiSi + m 2 J 2 S 2 + m3J353)5 = 
—mis — m^S — m^S = —m?S, 

SSJ = SJJ = -m^S, 

S{Sl + 51 + 532 ) = mi5Ji5i + m2SJ2S2 + m^SJ^S^ = miS^Ji + m2S^J2 + 
m352j3 = 52 J, 


SJS = {Si + 52 + 53 )(miJi + m 2 J 2 + rn 3 J 3 )S = (miJi 5 + m 2 J 2 S 3 + 

m3J3S2+miJiS3+m2J2S+m3J3Si+miJiS2+m2J2Si+m3J3S)S = J'52+ 
m2J2m3J3S + m3J3m2J2S + miJim3J3S + m3J3miJiS + mi Jim2J25 + 
m2J2miJiS = JS"^, 

S{Sl+Si+Sl) = {Si + S2+S3){Sl+Si+Sl) = 53 + 5| + 5f + 5i52 + 5i532 + 
S2Sl+S2Si+S3Sl+S3Si = -b‘^S-m3J3Si-m2J2Si-m3J3Sl-miJiSi- 
m2J2Si — miJiSi = —h'^S — msSfJs — m2SiJ2 — m3SiJ3 — mi52ji — 
m25|j2-mi53 Ji = - 62 ,S-( 5'2 + 52 _|_^ 2 ^ ^ 

-525 + (52 + 52 + Si)J + mi Ji52 + m 2 J 252 + m 3 J 352 = 

- b'^S + {Si + Si + Si)J + J 52 . 

Now we are ready to find the powers of (5 + 5). Using (21), we get 

{s+sf = 52 + 55 + 55+52 = s^+SJ+j s+sl+si+Si - SJ+j s = 
52 + 5J + 2JS - SJ + si + 5| + Si. 

Multiplying the result by 5 + 5 and applying again (21), we find 
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{S+Sf = {S+S)[S‘^ + SJ+2JS-SJ+Sf+S^+Si] = S^+S^J+2SJS- 
SSJ + 5(5f + Si + Si) + 552 + SSJ + 2SJS - S^J + 5(5^2 + 51 + 5|) = 
-525 + 52 J - 2m25 + m‘^S + S^J + JS‘^ + JSJ + 2SJS - [(5? + 5| + 5|) - 
SJ+JS]J+[-b^S-{Si+Sl+Si)J+JS‘^] = -(62+m2)5+252^+^52 + 

jsj+ 2 js^-isi+sl+si)j+sj‘^-jsj-b‘^s-{sf+sl+si)j+js‘^ = 

-ib‘^ + m2)5 + 252 J + 4^52 - 2{Sf + Si + Si)J - m^S - b'^S = 

- (62 + m2)(5 + 5) + 252 J + 4^52 - 2(5i2 + S^ + Si)J. 
Continue the process. 

(5 + 5)4 ^ (5 + 5)j_(52^^2)(5^5)^252j + 4^52-2(5i2 + 5i + 5|)^] = 
-(62 + m2)(5 + 5)2 + 2S^J + 45J52 - 25(5f + 5| + Si)J + 2SS‘^J + 
45^52 - 25(5? + 5| + Si)J = -(62 + m2)(5 + 5)2 - 2b'^SJ - Arr?S‘^ - 
2S‘^JJ + 2JS‘^J + 4^53 - 2[-625 - (5f + 5| + Sl)J + JS‘^]J = -(62 + 
m2) (5 + 5)2 - 2h‘^SJ - 4m252 + 2m'^S‘^ + 2JS‘^J- Ab'^JS + 2b'^SJ + 2(5? + 
5| + 5|) J2 _ 2JS‘^J = -(62 + m2)(5 + 5)2 - 2b‘^SJ - 2m‘^S^ - Ab^JS + 
2b^SJ - 2m^{Si + 5| + 5|) = -(62 + m2)(5 + 5)2 - 2m2[52 + SJ + 2JS - 
SJ + (5i2 + 51 + 5|)] + (2m2 - 262)(5^ + 2JS - SJ) = 

- (62 + 3m2)(5 + 5)2 + (2m2 - 262)(5^ + 2JS - SJ). 

Finally, 

(5 + 5)5 = (5 + 5)[-(62 + 3m2)(5 + 5)2 + (2m2 - 2b‘^){SJ + 2JS - SJ)] = 
-(62 + 3m2)(5 + 5)3 + (2m2 - 262)[52 J + 2SJS - SSJ + SSJ + 2SJS - 
52 J] = _(62+3^2j(^^^)3^(2m2-262)[52j-2m25-5^2+^5^+2^52- 
(512+51+51-5^+^5)^] = -(62+3m2)(5+5)3 + (2m2-262)[52^-m25+ 
JSJ+2JS‘^-{Sl+Sl+Si)J+Sj‘^-JSJ] = -(62+3m2)(5+5)3 + (m2- 
62)[ 252 ^ + 4^52 - 2(5i2 + 51 + Si)J - 2m^S - 2m‘^S] = -(62 + 3m2X5 + 
5)3 + (m2 - 62) [ 252 ^ + 4^52 - 2(5? + 5| + 5|)^ - (m2 + 62)(5 + 5) + (62 - 
m2)(5 + 5)] = —(62 + 3m2)(5 + 5)3 + (m2 —62)[(5 + 5)3 + (62—m2)(5 + 5)] = 
- 2(62 + m2)(5 + 5)3 - (62 - m2)2(5 + 5) 
which completes the proof. 


Thus, 

(5 + 5)5 = Lin (( 5 + 5 ) 3 , 5 + 5 ) ( 22 ) 

On the other hand, setting for brevity Rxy = R, from (19) we derive 


S + S=-R-2J = Un{R,J) (23) 

c 

Since (5 + 5) and J commute, (19) implies the commutation of R and J. 
Keeping this and j'^ = —m?E, from (23) we derive 

(5 + 5)3 = ^ i?3 + Lin {JR\ R, J) (24) 

(5 + 5)5= + Un{JR\R^,JR^,R,J) (25) 
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From (22), (23) and (24) 


{S + Sf = Lm Un{R^, JR^,R, J), Lin {R,J) =Lm{R^,JR^,R,J). 


So, finally from (25) 


R^ = Lin {JR^, R^, JR^, R, J). 


It is easy to trace that the coefficients of all linear combinations are 
polynomials in -,h, m. To complete the proof we should remark that 

^ Lin {JR^, R^, JR^,R, J)R = Lin (d^i^^ R^, JR^, R^, JR) = 
Lin [jLin {JR^, R^, JR‘^, R, J),R^, JR^, R^, JR 
Lin{R*,JR^,R^,JR,E) 
which allows to find all powers of R inductively. 
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